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Abstract 

Lepton flavor violation in r and /x processes is studied in the littlest Higgs model with T parity. 
We consider various asymmetries defined in polarized r and /i decays. Correlations among branch- 
ing ratios and asymmetries are shown in the following lepton flavor violation processes: ^'^ — )• e'^^, 
/i+ e+e+e", ^' A e~ A {A = Al, Ti, Au and Pb), r+ /i+7, r+ r+ ^+e+e", 

T+ — > /x+P (P = vr", r] and r?'), r+ — ;> /i+F {V = p^, uj and 0), r+ — ;> e+7, r+ e+e+e", 
r"*" — )• e"'"/x+/i~, r+ — )• e"^P, r"*" — )• e"''^, r+ — )• /u^^"^e~ and r+ — )• e"*"e"*'/u~. It is shown that large 
parity asymmetries and time-reversal asymmetries are allowed in /i+ — )• e'^e'^e" . For r lepton fla- 
vor violation processes, sizable asymmetries are possible reflecting characteristic chirality structure 
of lepton flavor violating interactions in this model. 
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I. INTRODUCTION 



Although the standard model (SM) is successful in describing almost all experimental 
results of high energy physics in terms of gauge theory, the dynamics responsible for the 
elerctroweak symmetry breaking is unknown. In the minimal version of the Higgs sector, 
one Higgs doublet field is introduced as an effective description of the dynamics below some 
cutoff scale A . Constraints from electroweak precision measurements indicate that A is at 
least beyond 3 TeV. On the other hand, this scale already introduces a fine-tuning problem 
of the Higgs mass term at the 0(1) % level. This is called a little hierarchy problem 

Little Higgs models ^| were proposed as a solution to the little hierarchy problem. In 
these models, the Higgs doublet field appears as pseudo Nambu-Goldstone (NG) bosons 
of new strong dynamics at the cutoff scale. A remarkable property is that the one-loop 
quadratic divergence to the renormalization of the Higgs mass term is cancelled by a proper 
choice of global and local symmetries, so that the cutoff scale can be pushed to O(10)TeV 
without a severe fine-tuning. A simple case is called the littlest Higgs model jst where the 
Higgs field is realized as a part of the NG bosons associated with the global symmetry 
breaking of SU(5) to S0(5) and the gauge symmetry is [SU(2) xU(l)]^. Then the quadratic 
divergence of the Higgs mass term in the SM is cancelled by the extra gauge bosons and the 
heavy partner of the top quark. Subsequent studies j4|, however, showed that these heavy 
partners need to be heavier than several TeV to satisfy severe constraints imposed by precise 
electroweak measurements, hence reintroducing some degree of the fine-tuning problem. In 
the littlest Higgs model with T parity (LHT) [5\ , the model is extended to have a Z2 parity 
so that the heavy gauge bosons assigned to be T-odd particles do not directly couple with 
a pair of the SM fermions, and the phenomenological constraints are somewhat relaxed. 

Flavor physics provides an interesting possibility to explore the LHT. In order to assign 
T parity, we need to introduce extra fermions (T-odd fermions) which are left-handed SU(2) 
doublets. Then, the flavor transition can arise in the vertices of heavy gauge bosons, T- 
odd fermions and ordinary SM fermions. There are two 3x3 unitary matrices describing 
these flavor transitions associated with quark and lepton sectors besides the ordinary fla- 
vor mixing matrices in the quark and lepton sectors i.e. the Cabibbo-Kobayashi-Maskawa 
(CKM) U and Pontecorvo-Maki-Nakagawa-Sakata (PMNS) j3] matrices. These matrices are 
new sources of flavor changing neutral current (FCNC) processes in the quark sector and 
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lepton flavor violating (LFV) processes in the charged lepton sector. In the early literature 
on this subject various observable quantities of FCNC and LFV processes were calculated in 
the LHT and showed that the effects of T-odd partners can be sizable and the correlations 
among various observable quantities can be different from other new physics model such as 
supersymmetric (SUSY) models was pointed out later that calculations missed 

a type of diagrams and the logarithmic dependences on the cutoff scale in the FCNC and 
LFV amphtudes disappear thanks to new contributions Q , thus reducing theoretical 
ambiguity associated with the physics at the ultraviolet cutoff scale. Branching ratios of 



various FCNC and LFV processes were reevaluated including the new contributions |15l. |16|. 
The qualitative feature turned out to be similar to the previous calculation though there 
were sizable changes at the quantitative level. 

In this paper, we present the results of further studies in tau and muon LFV processes. 
In addition to branching fractions, we also study observables deflned with the help of po- 



larizations of the initial muon and tau lepton [17|, ll8| . Polarized muon experiments can be 
done using surface muon beams that are 100% polarized in the opposite direction of the fi~^ 
momentum. In fact, a /i^ — > e^7 experiment with initial muon polarization is under consid- 



eration 



19| . In the tau pair production at e colliders, the polarization information can 



CERN 

cays 



be obtained by taking an angular correlation with tau decays on the opposite side [l8|. The 
study of polarized tau decays is also being considered for the r+ — )• ji'^ji^ji' process at the 
arge Hadron Collider (LHC) where we can utilize the tau polarization from W de- 
20|. In the following, we consider the processes n'^ — )■ e"^7, /i"^ — )■ e'^e'^e~ , fi^A — i- e~A 
{A = Al, Ti, Au and Pb), r+ /i+7, r+ r+ /i+e+e", r+ /i+P (P = 7r°, rj 

and T]'), r+ /x+V (V = u and 0), r+ e+7, r+ e+e+e~, r+ e+/x+/i", r+ e+P, 
r"*" — 7- e^V, r"*" — )■ fi'^fi^e^ and, r+ — )■ e^e^fi~ . We deflne a parity asymmetry in two-body 
decays and two parity and one time-reversal asymmetries in three-body decays. In addition, 
we can deflne forward-backward asymmetry and forward-backward-angular asjTumetries in 
the cases of r+ — )■ /i+e+e^ and r+ — )■ e~^fi'^fi~. We show that the parity asjTiimetries of 
two-body decays reflect the characteristic chirality structure of the LFV interactions. For 
three-body decays, we flnd that there are useful relations among various asymmetries. We 
calculate the rates of /i-e conversions for different muonic atoms and show that the ratios 
of the conversion rates can vary within 1 order of magnitude over most of the LHT model 
parameter space. These features, as well as correlations of various branching ratios, can be 
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useful in discriminating different new physics models. 

This paper is organized as follows: Section II is a brief review of the LHT. LFV processes 
are classified and various asymmetries are defined in Sec. III. Numerical results on the 
various observable quantities are shown in Sec. IV. Section V is the conclusion. Appendix A 
shows the general formulae of branching ratios and asymmetries for the processes studied. 
Functions for the Wilson coefficients are given in Appendix B. Useful formulae to perform 
the consistency test of the LHT are presented in Appendix C. 



II. THE LITTLEST HIGGS MODEL WITH T PARITY 



In this section we review the LHT in order to fix the notations we use in this paper. We 
use the Lagrangian given in Ref. Q with corrections discussed in Refs. QQ- 



A. Gauge and Higgs sectors 



Gauge and Higgs sectors of the littlest Higgs model are described as a nonlinear o model 
with the spontaneous global symmetry breaking from SU(5) to SO (5) with scalar fields, S, 
which is transformed as 15 representation of the SU(5). The vacuum expectation value of 
the scalar fields, Sq, which breaks the global SU(5) symmetry is 

/O 1 o\ 

1 

So = 1 . (1) 

1 
\0 1 / 

Generators of the unbroken S0(5) symmetry, T'^, and the broken ones of SU(5)/SO(5), X°- ^ 
satisfy the following relations. 







a\T 



a\T 



(2) 



The Lagrangian of the NG bosons and gauge fields is given by 

•^NG-gaugc 



£tr[(D^s)tD,s] + (Atr^rw.^;^ - \BrB,,^ , 

i=l,2 ^ ^ 



(3) 



4 



where / is the decay constant of the nonhnear a model and the NG boson field is written 
as follows, 



(4) 
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(5) 



Fourteen NG bosons are denoted by h, 7r^'°, r], 0^'°, 0^, and 0^^. As explained later, 
H = {-in+/y/2, {v + h + m^)/2f is identified as the SM Higgs doublet. V is the electroweak 



symmetry breaking vacuum expectation value, v = 246 GeV. and B, 



1, 2) are 



gauge fields of the SU(2)jXU(l)i, which are subgroups of the SU(5). The generators of the 
gauged subgroups are 





^0 ooy 

'^O ^ 



n=-diag(3,3,-2,-2,-2). 



(6) 



Q^2 





-a"^ 



r2=-diag(2,2,2,-3,-3). 



(7) 



where a" (a =1, 2, 3) are the Pauli matrices. The covariant derivative of E is defined as 



(8) 



i=i,2 



where g and g' are gauge coupling constants of SU(2)LxU(l)y. 



The Lagrangian 
is called T parity [5| 



3]) has a discrete symmetry under the following Z2 transformation, which 



n^-nnn, q = diag (1,1,-1,1,1). 



(9) 
(10) 
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The T parity eigenstates of the gauge bosons are 



V2 



W 



±,3At 



H 



^2 



(T-even), 
(T-odd), 



(11) 
(12) 



V2 ' ^ V2 
where Wf^ = (VT/^ t iWj^)/V2 (j = 1,2). As for the NG boson fields, we can show from 
Eqs. ([5]) and ffTOl) that h, 7r^'° are T-even particles and the others are T-odd particles. W^'^'^ 
and receive mass of the order of / from the gauge symmetry breaking [SU(2) xU(l)]^ — >■ 
SU(2)iXU(l)y by absorbing co^'^ and r]. After the electroweak symmetry breaking, 7r^'° 
are absorbed by the SM gauge fields, W^'^'^ and B'^. Mass eigenstates of the neutral gauge 
bosons are defined by 



wl" cose 



w 



B'i sin Ow, 



sin Our + B'^ cos 9 



(13) 

= W^j^ cos Qh - B^ sin Oh, = W^'' sin Oh + B>^ cos Oh, (14) 

where and are the SM Z boson and photon. 6w is the weak mixing angle determined 



as sin^i^ = g' / \/ ^ g'"^. The mixing angle of the T-odd gauge bosons is given by 

tan2gH = - . , J^'f:': (15) 



where = cos(f /(\/2/)) and = sm{v /{\/2f)). Expanding in terms of j p , we obtain 



sin Q 



H 



99 



y2 4/2 



+ 0(-) 



Gauge boson masses are given as 
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(16) 
(17) 

(18a) 
(18b) 
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1 _ _ 2 2 



5f^ 

4P+ lyi^ 



4 



^/2 



(18d) 



where sh = sin 9h and ch = cos 9h- We consider terms up to 0(f^//^) so that we neglect 
the difference between mw„ and m^^ in our analysis. 



B. Fermion sector 



The fermion sector of the LHT consists of three families of quark and lepton fields qi^2, 
^1,2, un, dji, uji, eji, qhr and £hb. for each generation, and a set of the "top-partner" fermions, 
t'^ 2 and t[j^2R- Following the procedure given in Ref. |l3|, we assign gauge charges of the 
fermion fields as summarized in Table HI The kinetic terms of left-handed fermions are given 
by 

+ ^^27^ [id, + V2gQfW^^ - V2g'B,^Y,^' - V2g'B2,Y,^') , (19) 
where these fermions are introduced as incomplete multiplets of the SU(5), 



V 













(20) 



and Pi^ = Under the T parity, they transform 



as 



Thus the SM leptons and the T-odd leptons are given as the T parity eigenstates, 

?i = — (T-even), 



V2 
V2 



(T-odd 



(21) 

(22) 
(23) 



respectively. The right-handed SM fermions are introduced in the same manner as in the 
SM. 

The right-handed T-odd leptons are introduced as a nonlinear representation of the SU(5), 
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(24) 
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TABLE I: Quantum numbers of fermion fields. Yi^2 are charges of U(l)i,2- Generation indices are 
suppressed. For qhr and Ihr, see the text. 



where upper three components, ipR and xRj heavy enough to be neglected in low energy 
dynamics, changes the sign with the T parity. In order to construct the gauge invariant 
kinetic term, we define as 



(25) 



which transform linearly under the SU(5) transformation. Thus the kinetic term is given by 



2^ - n I y- -A* • -J ■ ■ Ifj, 

+ (T parity conjuation). 



(26) 



where Y^" = ^diag(3, 3, -2, -2, -2) and Y^" = ^diag(2, 2, 2, -3, -3). The quark sector 
is constructed in the same procedure. For the quark sector, = ^diag(19, 19, 4, 4, 4) and 
F2*« = ^diag(16,16,16,l,l). 
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C. The Yukawa sector 



The Yukawa couplings of the T-odd fermions are 



^T-odd-Yukawa = -/^'Vl^S^ + n^O^^^^'n + H.C. 



with generation indices, i and j. Thus their masses are given by 



(27) 



m, 



1 + cJ 



m 



1 + 



(28) 
(29) 



where i is the generation index and are eigenvalues of the Yukawa coupling matrix, 
K*-'. We neglect 0{v'^/ p) differences between nf^jf^ and rn^fju calculation of LFV 

processes because these differences only contribute to the higher order corrections in the 
t>^/ p expansion. For the SM charged leptons, the gauge invariant Yukawa coupling terms 
are written as 



T-even- Yukawa 



feabe^yz {^^').{J:)ay{^)bz-{^f'^o)x{^)ay{^)bz 6^^ + H.C, (30) 



where a, b = (1,2) and x, y, z = (3,4,5) are SU(5) indices and i and j are generation 
indices. The left-handed SM lepton doublets are embedded in SU(5) multiplets \E'f * and 
\l>f * as 





V " / 



Xi 



( ^ 



(31) 



where X and X are SU(2)j singlet scalar fields whose U(l) gauge charges are (Fi, Y2) = 
(jq, —jq) and {—jq, jq), respectively. Following Ref. [2l|, we use (£33)^^/'^ and its T parity 
conjugate as X and X, respectively. 



D. New mixing matrices 

As explained in Ref. js], the mass matrices of the SM fermions and the T-odd fermions 
are not simultaneously diagonalized in general. Consequently, new mixing matrices are 
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introduced in the interaction among the T-odd gauge bosons and fermions. We extract the 
interaction terms of the T-odd W boson and leptons from Eq. (1191) as 



Acft D^{Vm)^,ly'HJ,PLi'W^'' + ^{VH,\ieHl,PLi^'W-^ + H.c. 



(32) 



where Vm and Vhv are new mixing matrices. These two mixing matrices are related to each 
other by 

VpMNS = V^^Vhu, (33) 

where Vpmns is the PMNS matrix because the mass matrices of in and uh are deter- 
mined by the same Yukawa couphng matrix k^^ . Vhi and Vhu also appear in the interaction 
of the neutral T-odd gauge bosons and leptons. We parameterize Vh£ as 



V, 
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5!2 
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(34) 



V 



where sfj = sinOfj and cf^ = cosO^j. We take O^j and 6jj as input parameters. Notice that 
there remain three unremovable phases in Vhi after the phase convention of the PMNS 



matrix is fixed, as pointed out in Ref. 



10|. For the quark sector, we define new mixing 



matrices Vnd and Vhu in the same way as Vne and Vhu- We take Vnd as a matrix parametrized 
ike Eq. f lM|) and the other mixing matrix, Vhu, is expressed by Vna and the CKM matrix 
6| in the following way: 



Vhu — yndVcxM- 



(35) 



III. LEPTON FLAVOR VIOLATION IN THE LITTLEST HIGGS MODEL WITH 
T PARITY 



In this section, we present 



In contrast to previous works 11 



'ormulae for observables in the LEV processes in the LHT. 



12 



14 



16| . we consider various angular distributions in 



polarized /i and r decays in addition to the branching ratios and the /i — e conversion rates 
in various muonic atoms. Studying the angular distributions is useful to extract information 
about the chirality structure of the low energy effective Lagrangian. 
The LEV processes that we consider in this paper are the following. 
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Radiative two-body decays: /i"*" — )■ e"'"7, r"*" — )■ /i+7 and r"*" — )■ e"'"7. 

Leptonic three-body decays of /i and r. We classify these decay modes into three types 
according to the lepton flavor combination in the final state. 

— Type I: /i+ — )■ e~^e~^e~ , r+ — )■ /i+yU+yU^ and r+ — )■ e+e+e^. All the three leptons 
in the final state have the same flavor. 

— Type II: r+ — t- fi~^e~^e~ and r+ — )■ e'^fi'^fi~. In these processes, the final state 
consists of a pair of same flavor leptons with opposite charges and another lepton 
with different flavor. There are no identical particles in the final state. 

— Type III: r+ — )■ fi^fi'^e'^ and r+ — > e^e^fi^. Lepton flavor changes by more than 
one. 

Semileptonic two-body decay of r. We consider r+ fi^P and r+ — )■ fi^V modes, 
where P and V stand for a pseudoscalar meson (vr^, 77 or r^') and a vector meson 
w or 0), respectively, as well as corresponding r"*" — e"^ modes. We do not consider 
r+ — i- /i+(e+)-ft'*-*'' modes, where both quark and lepton flavors are violated. 

fi — e conversion in a muonic atom. In these processes, fx^A^ e^A, we discuss 
aluminum (Al), titanium (Ti), gold (Au) and lead (Pb) as target atom. A, because 



22| 



the conversion rates depend on target nuclei 

In the LHT, LFV processes are induced by loop diagrams with the T-odd gauge bosons 
and leptons in the internal lines. Contributions of the T-even particles are negligible because 
they are suppressed by neutrino masses. Since only the left-handed SM (T-even) leptons 
couple to the T-odd gauge bosons, the low energy effective Lagrangian for LFV has a re- 
stricted chirality structure: the lepton flavor changing occurs in the left-handed lepton sector 
only. In the following subsections, we present the formulae for the effective Lagrangians, 
branching ratios, and various asymmetries, which can be defined based on angular distribu- 
tions. Except for the fi — e conversion, we explicitly show formulae for t to fi transitions. 
Those for /x to e and r to e are obtained by appropriate replacements of flavor indices. For 
simplicity, we neglect lepton masses in the final states. 

In the following, formulae of the low energy effective Lagrangians and the observable 
quantities are shown based on the LFV interactions relevant to the LHT. Those in the case 
of the most general low energy effective Lagrangians are given in Appendix |X1 
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A. — )• e''~7, — ?• ^'''7 and, — t- e^7 

The effective Lagrangian for r"*" — )■ yU"'"7 in the LHT consists of the following dipole 
moment type operator. 



A = - ^ [mrA'i'^fna^''^iLF^, + H.c] , (36) 



where 



AT = (4^^ Y.^VHdUVm)i2 [ncm{xi) + \nnm{xi) + ^iV^./(x;)) , (37) 



(47r)2 8/2^ 

with 



m^^ mj^^ tan^ 

The functions Ncm and A^atm are given in Appendix [Bl The coefficient of the operator with 
the opposite chirality, fia^^ HfiF^i,, is suppressed by m^/mr due to the chirality properties 
in the interaction of the T-odd gauge bosons. 

For a r+ — 7- yU'*'7 decay of a polarized r"^, we define a decay angle 9 as the angle between 
the spin of r+ and the momentum of in the rest frame of . The configuration of spin 
and momenta are shown in Fig. [H 

The differential branching ratio with respect to cos 9 is given by 

-^Br(r+ ^ /x+7)lht = r.^^ \A'^r\' (1 " ^os^). (39) 
a cos 9 n ' ' 

We define the polarization asymmetry, which characterizes the chirality structure of the 

effective interaction in fl36l) . as 

where the integrated branching ratio is 



Br(r+ ^ /i+7) = ^ 



^rfBr(r+ ^ /i+7) 

cos^ — 41 

ct cosy 



is a parity odd quantity since cos 9 is proportional to the inner product s ■ p, where s 
and p are spin of initial and spatial momentum of final /i"*", respectively. By substituting 
fl39|) . we obtain 

Br(r+ ^ /.+7)lht = r.^^ , (42) 

A^(r+ ^ /i+7)LHT = (43) 
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The parameter independent result of is a consequence of the special chirality structure 
of the LHT. 




FIG. 1: Configuration of — t- /i^7 decay. The dashed arrow denotes the direction of the r"*" 
spin. The solid and the wavy arrows show the directions of momenta of the fx~^ and the photon, 
respectively. 6 is dehned as the angle between the directions of the r"*" spin and the iJ,~^ momentum. 

B. — )• e^e'^e", r+ — )• and r+ — )• e+e^e~ 

The effective Lagrangian for type I leptonic three-body decay, r+ — )• in the 

LHT is given as 

AC 

+ ^i^"^(rL7'^/^L)(/xWij) + H.c] , (44) 
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where 



2 „,2 



I,LHT _ 9 'V V^/T^ X* /T^ X 



(47r)2 8/2 



I 



- sin^ 6'iyP^(xi) + ^-^ + sin^ -Pz(a;0 + ^Y^{VHi)*„a{VHi)m2B(^e){xh Xm 



(4vr)2 8/ 



The functions P-^, Pz and P(e) represent contributions of photon penguin, Z penguin and 
box diagrams, respectively, and are given in Appendix |Bl A^"^ is shown in (137|) . 

The spin and momentum configuration of a three-body decay of a polarized r"^ at rest is 
determined by two energy variables and two angular variables, as depicted in Fig. [2l The 
three spatial momenta of the particles in the final state are denoted as Pa, Pb, and Pc- We 
take z axis in the direction of Pa and the decay plane is identified as the z-x plane. The 
direction of x axis is chosen in such a way that the x component of pb has a positive value. 
We take the y axis along pa x Pb- The direction of the polarization vector of the initial 
particle is parametrized by the polar and the azimuthal angles 6 and 0, respectively. For the 
r+ — i- decay, we identify the momentum of /i~ as pa- Two /i"^'s are distinguished 

by their energies. The momentum of the /x+ with larger (smaller) energy is identified as pb 

m. 

As for the energy variables, we define 

Xb,c = , (47) 

rriT 

where Eb^c are energies of /i^'s with momenta pb^c- Eb is larger than Ec by definition. Possible 
ranges of Xb,c in an approximation of massless final particles are 

^ < X6 < 1, 1 - Xb < Xc < Xb- (48) 

However, the r+ — )■ fi^fi^fi^ decay amplitude is singular for Xb ^ 1, where the invariant 
mass of the pair of yU"*" (with p^) and fi~ is vanishing within the approximation of neglecting 
the muon mass. The singularity comes from the pole of the photon propagator in the 
contribution from the dipole moment operator. We introduce a cutoff parameter < 6 1 
and take the integration interval of 

^<Xb<l-5, (49) 
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in order to avoid this collinear singularity. 




FIG. 2: The spin and momentum configuration of a three-body decay in the rest frame of the 
initial lepton. The decay plane is identified as the z-x plane where the z axis is taken to be the 
direction of particle a's momentum, the x axis is defined so that the x component of particle 6's 
momentum is positive. The direction of the initial lepton polarization vector s is parametrized by 
the polar and the azimuthal angles 6 and (f), respectively. 

The (5-dependent branching ratio is calculated by integration over Xb,c, cos 6 and as 

Br(5)= / dcos9 dxf, ^^a;,— — — — — . (50) 



1 Jo J 1/2 Ji-Xb dxbdxcdcj) dcosO 

The general formulae of differential branching ratios are shown in Appendix |Al We also 
define the following angular asymmetries: 

^^('^) = ^5"77T / ~/ ]dcos9 d(t) I dxb dxc-^ — -j—: -, (51a) 

Br(()) \Jo J^iJ Jo A/2 Ji^x, dxbdXcdcpdcosO 

1 ( /"'^/^ /"Stt/sX rl-& rxb d'^Bl 

^^('^) = W7J\ dcosO [ - ]d(t) dxb dx^-—^ — -j—. (51b) 

Br((3) y_i \J-k/2 Jn/2 J Ji/2 Ji~x, dxbdXcd(pdcos9 

Ay{5) = / dcos6 { - ] d(j) dxb dxc-: — i — j—: (51c) 



Br((5) J_i \Jq J-ttJ J 1/2 Ji-x^ dxbdxcdcj) dcosO 

The branching ratio and asymmetries in the LHT are expressed as 

Br(r+ ^ /iV/i")LHT(5) =Br(r+ ^ z/.e+z/e)5i'^"^(5), (52) 
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C^'rM4W + C^t^^AW), (54a) 



^^""^('^) = ^^i.LHT(^) (-Cj^;-"^Ci((^) + C'^^C,{5) - C™C3(5) + CS^^C4(5)), (54b) 
^'^^l'^) = 2^i,LlT(^) (-g''5"^^3(^) + CT^CS))- (54c) 

The functions Ai(5), Bi{6), and Cj((5) are given in Appendix lA 2[ Cj^^^""", etc. are written 
in terms of the Wilson coefficients in (HH) as 



^I,LHT _ 1^ /iLHT|2 ^I,LHT 



I,LHT 



^I.LHT 



I,LHT 

9u 



(55a) 

Cit"^ = ReMr^;;^"*], C^S"" = Re[eA--,i-^*], (55b) 

= imMr^^^n, = imMr^L^^^i. (55c) 

As explained in Appendix IA21 Az and are parity odd asymmetries and Ay is a time- 
reversal asymmetry. 

The time-reversal asymmetry Ay is induced by the CP violation in the effective La- 
grangian, since we are considering a theory with CPT invariance. In fact, Cj^^""" and Cj^^"^ 
are proportional to the Jarlskog invariant [23'] defined for Vne'- 



CjIjI 0^ JiVm) = Im [iVm)niVmy22iVm)i2{Vm)2i] 

Therefore, a three-generation mixing is necessary for nonvanishing Ay. Notice that two 
of three complex phases in Vhi are removable when T-even neutrino masses are neglected. 
That is why J{Vhi) depends on only one combination of the phases. 
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C. r"*" — )• fj,~^e~^e and r"*" — )• e^n^^fi 

The effective Lagrangian for a type II trilepton decay r+ — )• yU+e+e~ has a similar struc- 
ture to fHIjl : 

+ ^?";'^"^(rL7>L)(eij7^e«) + H.c] , (57) 

where 

^S"""^ = (1^^ (E(^^^)*3(^h.).2 [- sin^ OwP.ixi) + + sin^ ew)Pz{xi) 



l,m 

+ {Vm)UVHe)ii{VmyrniiVmU2\Bie){xi,x^)y (58) 

and g^ll^^^ = 9^1^^'^ ■ ^)^^ is the same as (137|) . 

Kinematical variables are defined in a similar way to the type I case. For the momentum 
assignment, we identify particles a, b and c in Fig. [2] as /i"*", and e~, respectively ^. A 
difference from the type I case is that both Xb > Xc and Xb < Xc are possible in type II modes 
because there are no identical particles in the final state. Therefore, the range of Xb and Xc 
are 

6 <Xb<l, l + 6-Xb<Xc<l, (59) 
where the cutoff parameter < 5 ^ 1 is introduced to avoid the collinear singularity in 

Xb + Xc ^ 1. 

Formulae for the branching ratio and the angular asymmetries Az^x,y in ~^ fi^e^e~ 
decay in the LHT are obtained as 

Br(r+ ^ /x+e+e-)LHT(5) =Br(r+ ^ z7.e+z/e)i?"'^"'"(^), (60) 
+ (C5^^^ + C™)D4(5), (61) 



^ The convention of particle assignment is different from that in Ref . 18 1 
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-i(C™ + C?3--)D,(5)), (62a) 

^"■^""(^) = ((-C^S'"" + C^i^f + (C^S'"" - C'r')Em , (62b) 

^"''"'(^) = ^^ilW ( - (^S"' - ^ir")^3(5)) , (62c) 
where 

^II.LHT _ I ^LHT|2 ^II,LHT _ II,LHT ^ ^II.LHT _ II.LHT ^ /^o \ 

^Rl — |<^^i?. I ' '^L3 — dhr ' *-"L4 " 9 LI ' I^Ddaj 

= MeA-^-gir\ = Re[eA--,£--], (63b) 

CS-^ = ImMr.S^^^I, C^ir" = ImMr^i^;^""*]- (63c) 

The functions Di{6), Ei{6) and -Fi(5) are given in Appendix lA 31 

In the type II decay, another class of asymmetries can be defined. We define the energy 
asymmetry in r+ — )■ yU+e"*'e~ as the asymmetry between the partial widths with i?e+ > E^- 
and Ee+ < E^- in the rest frame of initial r"*": 

^""""^^Tif "/ ) '^''''^'''l^- ^^^^ 

'-'^ \Jxf,>Xc Jxf,<Xc/ UXfoUXc 

This asymmetry is also called as the forward-backward asymmetry, because Xb > Xc {x^ < Xc) 
corresponds to -pe- > (p^+ -pe- < 0) in the rest frame of the e"*" e~ pair. Furthermore, 
we define the following asymmetries combining asymmetric integrations over [xb, Xc) and 
(cos 6', (p): 

AzFB = „ , ( [ — [ \ dcosO [ d(j) ( [ — [ 1 dxhdx. 



Br(5) VJo J-iJ Jo \Jx,>x, Jx,<xJ dxbdxcd(f)dcose 

(65a) 

^^PB = i^K^j dcoself -[ ]d(j)([ -[ ) dxbdxc ^ ^, 

Br(5) J_i A/2 / \Jx,>x, Jx,<xJ dxbdxcd(j)dcose 



(65b) 



(65c) 



^^^^ = ^5"7xY / ^cos^ / - / ]d(p{ - \dxbdxc-r 

Br((5) J_i \Jo J_^J VA,,>xe Jx^<xJ dxbdxcd(t)dcose 
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In the LHT, we obtain 

"^FB (^) = ^II.LHT I ~ I^^LZ -Cla )D2{d) + -{Cj{ - C jDiid)], 

(66a) 

^">T(^) = - C'lDD.iS) - liC-r- - C-i^nDAiS)^ , (66b) 

+ ^(CS^«^ + C™)E2(5)), (66c) 
= ^^inr ( - ^(Cj7 + Cy,)E,{5)^ , (66d) 
where the functions Gi{6) and Hi{6) are given in Appendix lA 31 



D. r+ —7- fi'^ fj.'^e and r+ — t- e^e^/i 

The type III leptonic decay modes, such as r+ — )■ fi^fi^e^ , are those in which the lepton 
flavor changes by two. At the one-loop level, only the box diagrams contribute to these 
decay modes. The effective Lagrangian for r"*" — )■ fi~^fi~^e~ is 

An = [^7"/''"^(rL7>L)(ei7./^L) + H.c] , (67) 

where 

III LHT 9^ ^ 

^m.LHT = ^ ^{VHeyi3iVm)i2{VHiymi{Vm)m2B^e){xi,x„^). (68) 

^ ^ ■' l,m 

The kinematics of the type III decay is treated in a way similar to the type I case. 
For r"*" — )■ fi^fi^e^ mode, e~ and one of the /x^'s, which has a larger energy in the rest 
frame of r"*", is identified as particles a and b in Fig. [2l respectively. We define the angular 
asymmetries Az, Ax and Ay by fISlap . (ISlbp and fISlcp . respectively. The forward-backward 



type asymmetries are not defined because there are identical particles in the final state as 
in the type I case. 

The branching ratio in the LHT is written as 

J^f^'^y^^^^ =^^S'^""A3(^), (69) 
Br(r+ -> z/^e+z/e) ^ ^ ^ ^ 
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where 

(70) 



^III,LHT 
^L3 



III.LHT 

9li 



and ^3(5) is given in Appendix IA2[ Since the differential width is written in terms of one 
combination of the Wilson coefficient, C™'^^""", the angular asymmetries are determined 
independently of the input parameters: 

^ni,LHT ^ _ 1 ^ ^ni,LHT ^ ^in,LHT ^ (7^) 



E. — )• TT , r] , 1]' and r+ — )• e^TT,r],r]' 

The effective Lagrangian for semileptonic r — > decay processes are written as 
4Gf 



+ 5Z (^?LJ(J(^L7>L)(gL7/.gL) + (7i"(g)(rL7''/i^ (72) 



The Wilson coefficients are given as 
3m^) = (47r)2 8/2 5Z(^^^)«(^^^)'2 



X 



(5gSin^6'vK-P7(a;z) + (T^ - Qg sin^ 6'vK)-Pz(a;/) + ^{VHq)*mi{yHq)miB(^q){xi, Xm) I , 

(73) 

2 2 

^SS = ^^E(^^^)«(^^^)'2Q''Sin2^H^(^7(^0 -i^z(xO), (74) 

where and Qg are the weak isospin and the electromagnetic charge for quarks, respectively, 
is given in (ISTj) . 

For the two-body decay process r+ — i- fi^P, where P stands for a neutral pseudoscalar 
meson ir'^ , rj or rj' , the differential branching ratio is written as 

rfBr(r+ ^ /i+P)LHT ^ G|m^ / V lr^LHT|2 , 

3 Z = ^ 2" I'^RP \ (l + cost^), (75) 

where the decay angle 6 is defined in the same way as in the r+ fi'^'j case. The effective 
coupling G]^ for P = tt^ and P = rj are written in terms of the Wilson coefficients in (172|) 
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as 

^LHT _ /tt / LHT LHT LHT , „LHT\ frjf>\ 

~ 2v^ ^ ^''("^ ~ 9Lr{d) ~ 9lI{u) + 9Ll{d)) 5 \'^) 

^LHT _^ ( fv / LHT , LHT _ LHT _ „LHT N , fsf^LHT _ „LHtA 

'^Rt] — 2 \~y^^9Lr{u) ~^ yLr{d) 9 Ll{u) 9 Ll{d)) + J r)\9 Lr{d) 9Ll{d))j- {' ' ) 

and are the decay constants of the pseudoscalar mesons that are defined in Ap- 
pendix lA 5[ The effective couphng for r"*" — )■ fi'^rj' is obtained by replacing the subscript 
7] with 7]' in the expression of Gjirj. Formulae for the branching ratio and the polarization 
asymmetry Ap are derived from fl75|) as 

Br(r* ^ = (l - g)' |G«pI^ . P8) 

/1p(t+ ^ fi+P)LHT = i. (79) 

is purely determined by the chirality structure of the effective Lagrangian ( !72|) and 
independent of the values of the input parameters. 



F. — )• fi'^p, oj, <j) and r"*" — )• e"^/?, lo, (j) 

The semileptonic two-body decay r"*" — )■ where V denotes a neutral vector meson 

p°, cu, or are also described by the effective Lagrangian fl72|) . The branching ratio and the 
polarization asymmetry Ay-, which is defined similarly to A^ and Ap, are written as 

Br(r+ ^ /z+1^)lht = f 1 - 4 V ^v""^' (80) 

vr V mi 



By^^ — \Grav\'^ ( 2 + + I^'lvI^ — ^ Am^ — ~ ~ '^^^^^ ^^^^*lv\^ (^1) 



jLHT _ \r^^^ |2 / o I '"V 111/-- |2 '"r ^"^y 

1 



Ay(r+ ^ /i+\^)LHT = 25LHT 



(82) 
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The effective coupling constants for V = p'^, u and are 



^LHT 



/-yLHT 



yoLHT 



^LHT 



^LHT 
'-^RActi 



y/2mp 

fp'^^P / _LHT „LHT I „LHT „LHT \ 
2y^m ^ "^"^ yLl{d) ~^ yLr{u) yLr{d)J 



^^R ) 



/^"^^ ('^LHT , „LHT I ^LHT , ^LHT \ 
2y/2m ^ yLZ(d) + yLr(n) + 9Lr{d)) 

f<t> "^T- ^ /iLHT 

3 



^LHT 



2mr '^^^^^'^ ^ ^^^'^"^'^ ' 
where the decay constants fp,uj,ii) are defined in Appendix lA 51 



(83a) 
(83b) 
(83c) 
(83d) 
(83e) 
(83f) 



G. II — e conversion 



The effective Lagrangian for the coherent p—e conversion processes has the same structure 
as ([72]). 



had 



pv 



q=u,d 



The Wilson coefficients are obtained by replacing the lepton flavor indices 3(r) and 2(yu) 
with 2(/i) and 1(e) in Eqs. (1371). (I73|) and (17^1). 



We calculate the p — e conversion rates following the method given in Ref. 22|. The — e 
conversion branching ratio in the LHT is written as 



R(/i A^e v4)lht 



2GI 



'-'^capt 



-A^^D + 2(2,^HT + 2,i;HT + ^lht + ^lht )^(,) 



, r,/ LHT , LHT , O^LHT , o„LHT ^T.(n) 
+ "^{dniu) + 9Lr(u) + "^duid) + '^9Lr{d)) ^ 



r(u) 



(d) 



'■(d)' 



(85) 



where u 
in Ref. 



is the muon capture rate. y*^^-* and V'^'^'^ are the overlap integrals defined 
22| . For the reader's convenience we quote the values for Al, Ti, Au and Pb in 



Appendix lA 61 
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IV. NUMERICAL RESULTS 



In this section we show numerical results for the LFV observables given in the previous 
section. We have to specify 19 parameters in order to calculate LFV observables in the 
LHT: the decay constant /, six masses of the T-odd fermions and m'^f^q, six mixing 
angles and six phases in the mixing matrices Vne and Vnd- The value of / is constrained 
to / > 500 GeV by the electroweak precision measurements [24]. Throughout the analysis 
in this paper, we fix / as / = 500 GeV. Since all the LFV amplitudes are proportional 
to as shown in Sec. Illlt the branching ratios scale as f~^. We also assume that the 
T-odd quarks (except for the top partner) are degenerate in mass for simplicity. Under 
this assumption, T-odd particle loops do not induce additional contributions to the quark 
FCNC observables, and the mixing and the phase parameters in Vnd are irrelevant. We 
fix the T-odd quark mass as 500 GeV. The top-partner quark mass is irrelevant for the 
LFV studied here. There remain nine free parameters in the T-odd lepton sector. We vary 
these parameters independently within the ranges 100 GeV < < 1 TeV {i = 1, 2, 3), 
< e^,j < 271 [ii]) = (12), (23), (13)] and < ^ < 27r [(zj) = (12), (23), (13)]. For the 
type I and type II trilepton decays, we use following cutoff parameters: 

2 



3 X 10-^ ~ 3 (^^)^ e+e+e", 
0.04- 3 :r+^/x+/i+/x-, r+^e+/i+/i", (86) 

10-6 ^ 3 f^)^ : r+ ^ e+e+e", r+ ^ /i+e+e". 



For the type III modes we take 5 = since there are no singularity in the differential widths. 

We evaluate the observables explained in Sec. Illll for each model parameter set, and check 
if that set is allowed under current experimental constraints. At present, the upper limits of 
the branching ratios are available for various LFV processes. We use the values summarized 
in Table HI 



A. ^1 LFV 

The current experimental bounds of /i+ — )■ e+7 and /i+ — )■ e+e+e" are given by the MEGA 
25j and SINDRUM I |26j collaborations, respectively. Also, the SINDRUM II collaboration 



provides us with upper limits on the /i — e conversion rates for titanium [27], gold 
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Mode 




Upper limit Ref. 


Mode 








Upper limit Ref. 


— )> 


"7 


1.2 X 10-^^ 


[251 


r+ ^ 




'7 




4.4 X 10"*^ 


[301 


u+ — > 




1.0 X 10-12 


[26] 


r+ ^ 


1 


>^ 

1 


1 


3.2 X 10-8 


[31] 

[-■■ ■■- .J 


;U~Ti— )• 


e-Ti 


4.3 X 10-12 


[27] 


r+ ^ 




-e+ 


e- 


2.7 X 10-*^ 


[31] 




e-Au 0.7 X 10-12 


[28] 


r+ ^ 








1.1 X IQ-'^ 


[32] 




e-Pb 


4.6 X 10-11 


[29] 


r+ ^ 


/^^ 






6.5 X 10-8 


[33] 










r+ ^ 


/^^ 






1.3 X lO-'^ 


[33] 










r+ ^ 








2.6 X 10-8 


[34] 










r+ ^ 








8.9 X 10-8 


[35] 










r+ ^ 








1.3 X lO-'^ 


[35] 


Mode 




Upper limit Ref. 


Mode 








Upper limit Ref. 


T+ — ;> 




3.3 X 10-^ 


[30] 


r+ ^ 




V" 




2.3 X 10-8 


[31] 


r+ 




3.6 X 10-s 


[31] 


r+ ^ 


e^ 






2.0 X 10-8 


[31] 


T+ — ;> 


>+^- 


3.7 X 10-^ 


[36] 














T+ — )• 




8.0 X lO-'^ 


[33] 














r+ — > 




9.2 X 10-^ 


[33] 














r+ — ^ 




1.6 X lO-'^ 


[33] 














r+ e"* 


>° 


4.6 X 10-s 


[34] 














r+ — ;> 


'UJ 


1.1 X lO-'^ 


[37] 














T+ — ;> 




3.1 X 10-^ 


[34] 















and lead 



TABLE II: Experimental upper bounds of various LFV branching ratios. 
29| . The upper bounds given by these experiments are summarized in Table Ull 



At present, the MEG experiment is ongoing in search for /z"*" — )■ e"'"7, and an upper bound 



Br(/x"'" — )■ e"'"7) < 1.5 x IQ-n has been reported recently 38|]. For yU — e conversions, COMET 
and Mu2e experiments are in preparation jsQ]. 

We present correlations among the branching ratios of /i"*" — )■ e'*'7, e^e+e- and 

//-Al — 7- e-Al in Fig. [31 The color of each dot represents the value of Br(yU+ — t- e+7): 
black, red/gray and yellow/light-gray correspond to IQ-i^ < Br(/i+ — )■ 6^7) < 1.2 x IQ-n, 
IQ-i^ < Br(/i+ e+7) < 10"i2 and, Br(/i+ e+7) < IQ-i^, respectively. In all the 
scatter plots hereafter, we use the same color code. The current experimental upper limits 
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f -5DDQbV, m'^- sm^V, ImCkV ITbV, e^.^jiarbltrGir 



f -SDDQaV, m|^- SDCkV, ImCkV ITbV, A^.^jiBrbltmr 



b 



CD 



b i 

r- 

'b 

n 

15 



p 1 1 1 1 r 

in-'' ^a-'" id-'" iir'* 10"'= iir'^ lO"" 
Br{ -^e'y] 




ID"" in-'" in-'= in-'* id-'" in-'= id-" 



(a) 



(b) 



f -5DDQbV, m'^- sm^V, ImCkV ITbV, e^.^jiarbltrGir 



1 
< 



b 

HI 

'b 

Lfl 

'b 

n 

15 




— p — I — I — I — I — 

in-'= iQ-" ID-" iQ-" ID-" in-'^ ID-" 
Br( [i"^ e" ) 

(c) 

FIG. 3: Correlations among the branching ratios of /i"*" — )• e"'"7, /U"*" — )• e~^e~^e~ and ;U~A1 — J- e~Al 
in the LHT, where the decay constant / is taken as 500 GeV. The T-odd lepton masses are varied in 
the range from 100 GeV to 1 TeV. The mixing angles and the phases in the mixing matrix Vne are 
also varied in the whole range. The T-odd quark masses are fixed as 500 GeV. The horizontal and 
vertical lines are experimental upper bounds. The dashed lines show the branching ratios calculated 
with the dipole contributions only. The color of each dot represents the value of Br(^^ — t- e~^j): 
black, red/gray and yellow/light gray correspond to 10^"*^^ < Br(/i+ — t- e+7) < 1.2 x 10^^"*^, 
10"^^ < Br(/z+ e+7) < 10"^^ and Br(/x+ e+7) < 10"^^, respectively. 
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on Br(yU"'" — )■ e'^'y) and Br(/i"'" — )■ e"'"e"'"e~) are shown by horizontal and vertical lines. The 
dashed lines show the branching ratios calculated with only the contributions from the dipole 
moment type operator Jj^b.^^^ ^lF^u- In a certain class of models, such as SUSY, it is known 
that the LFV effect dominantly appears in the dipole moment operator. In such a model, 
the branching ratios are predicted to be on the dashed lines in the correlation plots. In 
contrast with those models, the correlations in the LHT show that the contributions of the 
dipole moment operator are less significant than those from other four-Fermi interaction 
terms in /i"*" — )■ e^e^e~ and /i^Al — > e^Al, as noted in Ref. 12(|. We can see that there 
is rather strong correlation between the branching ratios of /i^ — ?■ e"'"7 and yU^ — )■ e^e^e" 
whereas no correlation is observed in the branching ratios of ii^ — )> e^7 and the conversion 
rate for Al. 

We show the angular asymmetries of yU^ — )■ e^e^e^ in Fig. HI As explained in Sec. IIIIBI 
Az and Ax are parity odd asymmetries and Ay is a time-reversal asymmetry. We can see 
that Az is within the range from about —45% to +30%, while Ax is within the range from 
about —10% to +20%. There is a positive correlation between Az and Ax as shown in 
Fig. 111(c). Both asymmetries do not have correlations with the branching ratio. Possible 
value of the time-reversal asymmetry Ay is in the range from about —10% to +10%. 

In Fig. [5l the /i — e conversion rates for Al, Ti, Au and Pb are plotted as functions of the 
mass of the first generation T-odd lepton m\j^ for / = 500 GeV. Here, we fix the parameters 
in the T-odd lepton sector as m\j^ = m]j^ = 400 GeV, = 7r/500, = = and ^ = 0. 
T-odd quark masses are also fixed to 500 GeV, as in other scatter plots in this paper. The 
region > 500 GeV is excluded since the branching ratio for Au exceeds the experimental 
upper limit. At m}^^ = 400 GeV, all the three T-odd leptons are degenerate in mass, so 
that all the LFV amplitudes vanish. We can see that the conversion rates are suppressed 
also in a region between = 100 GeV and 200 GeV. In this region, cancellation among 
transition amplitudes occurs at different points for different nuclide. 

We show correlation between the fi — e conversion rates for Al vs Ti (Au, Pb) in Fig. [6l 
We can see the ratios of the conversion rates vary within 1 order of magnitude in most of 
the parameter space. We also notice that, in some cases, the conversion rate for Ti, Au and 
Pb can be close to the experimental bounds even if the rate for Al is suppressed. 
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f -5DDQbV, m'^- SIAV, ImCkV <m'„< ITbV, e^.^jiarbltrGir 





(c) (d) 

FIG. 4: Angular asymmetries of /i"*" — )• e~^e~^e~ as functions of the branching ratio in the LHT 
for the same parameter set as in Fig. [3l The correlation between Az and Ax is also shown in (c). 
The vertical solid lines in (a), (b), and (d) is the experimental upper limit of Br(//+ — )• e"'"e+e~). 
The color code is the same as in Fig. [3l 
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f = 500GeV, m'nq= 500GeV, lOOGeV < < 1TeV, m^^f= 400GeV, e',2=7i/500 




100 200 300 400 500 600 700 800 900 1000 

mi, [GeV] 

FIG. 5: jjL — e conversion rates for lead, gold, titanium and aluminum as functions of the first 
generation T-odd lepton mass ?ti^^ for / = 500 GeV. Other parameters in the T-odd lepton sector 
are fixed as m^^^ = m|^^ = 400 GeV, 6*^2 = vr/500, 6*^3 = 6*^3 = and 5^^ = 0. The T-odd quark 
masses are also fixed as 500 GeV. The horizontal lines are the experimental upper bounds for Pb, 
Au and Ti. 

B. r LFV 

The current upper bounds for r-LFV decays are listed in Table Ull These bounds are set 
by either the Belle and the Babar experiments. Improvements by 1 or 2 orders of magnitude 
are expected at future B-factories at KEK and in Italy [4^. 

In this subsection, we mainly present the results on the observables in r — |U decays 
discussed in Sec. Illli Quantities in r — )■ e decay modes behave similarly to corresponding 
ones in r — )■ /i modes. Correlations between the observable quantities in r — /i and r — )■ e 
modes are discussed in Sec. IIV CI 
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f - 5DDQbV, m'^- SIAV, 1 mCkV <m'„< ITbV, e^.^jiarbltrGir f - SDDQaV, m|^- SDCkV, 1 mCkV <m'„< ITbV, A^.^jiBrbltmr 




R( [1- Al -> e" Al ) R( [I- Al -> g" Al ) 

(a) (b) 



f -5DDQbV, m'^- SIAV, ImCkV <m'„< ITbV, e^.^jiarbltrGir 




R( [1- Al -> e" Al ) 

(c) 

FIG. 6: Correlations among the /i — e conversion rates for Ti, Au, or Pb vs Al. The horizontal 
lines are experimental bounds. The input parameter set and the color code of the plots are the 
same as those in Fig. [3l 

1. T ^ fi^ and trilepton decay modes 

In Fig. U\ we show correlations among the branching ratios of — )■ /U^7, — fi^fi^fi~, 
r+ —7- yU+e+e" and r+ — )■ fj^^fi^e^. We see that the branching ratio of r+ — )■ fi'^'j can be 
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f -5DDQbV, m'^- SIAV, ImCkV <m'„< ITbV, e^.^jiarbltrGir 




in-'" iQ-^' ID"" ID-'- iiT ID"' 



ID"'" in-'' in-" in-"" in^ in^ iir' 



(a) 



(b) 



f -5DDQbV, m'^- sm^V, ImCkV ITbV, e^.^jiarbltrGir 




in-'" IQ-" ID-" IQ-"" Iff^ Iff^ ID-' 
Br(T- ) 

(c) 

FIG. 7: Branching ratios of r"*" — )• r+ — )• ^"'"e^e" and r"*" — )• fi'^fi^e~ as functions of 

Br(r^ — ?■ /x^/i+yU^) for the same input parameter set as in Fig. [3j The horizontal and the vertical 
lines are the experimental bounds. The dashed lines show the branching ratios calculated with the 
dipole contributions only. The color code is the same as in Fig. [3l 
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as large as 10~^, which is close to current experimental upper limit. For r"*" — t- jX^ fi^ 
and r+ — )■ /i"'"e"'"e~, the possible maximal values of the branching ratios are about 1 order 
of magnitude below the corresponding experimental limits. The behavior of the correlation 
between Br(r+ — )■ yU^7) and Br(r+ — )■ fi^^^fi^) shown in Fig. [7](a) is similar to the fi e 
case given in Fig. [21(a): Br(r+ fi^^^fi^) is larger than the prediction in the dipole- 
dominant case (dashed line) by 1 or 2 orders of magnitude. Since we take all the masses 
and the mixing angles/phases in the T-odd lepton sector as free parameters, there is no 
direct correlation between the r — )■ and — )■ e transition amplitudes. Therefore, the 
correlation plots among r — t- /i processes do not change much even if the upper limit of 
Br(/i — i- 67) is lowered. As shown in Fig. [7](c), the branching ratio of r+ — )■ /i+/i+e^ is 
highly suppressed because this process involves the /i — e transition, whose magnitude is 
constrained by other yU — t- e processes discussed in the previous subsection. We see that the 
suppression of Br(r+ — t- fi^fi^e^) is stronger for a smaller value of Br(/i"'' — i- e'^'y). As for 
the other type III decay mode r+ — t- e+e+yU~, the correlation plot with Br(r+ — )■ e~^e~^e~) is 
almost the same as Fig. H^c). 

We show angular asymmetries of — i- fi^fi^^~ in Fig. [HI We find that the parity asym- 
metries are within the ranges —25% < Az < +25% and +5% < Ax ^ +15% irrespective 
of the branching ratio for Br(r+ — )■ yU+/i+/i~) > 10~^^. The time-reversal asymmetry Ay 
is very small for the same range of the branching ratio. Compared with Fig. [U we can see 
that the possible ranges of the asymmetries for r+ — )■ /i^/i^/i" are narrower than those for 
/i"*" — )■ e~^e~^e~ . This quantitative difference between the r+ — ji^ jj,^ ji^ and /i"^ — )■ e"'"e"'"e~ 
cases is caused by the difference in the allowed ranges of the branching ratios. In the 
general scan in the parameter space of the T-odd lepton sector, the decay amplitudes for 
/i"*" — )■ e^e^e~ , r+ — )■ yU+/i+yU~ and r+ — )■ e+e+e^ behave in the same way. Consequently 
the distribution patterns of the asymmetries in the scatter plots such as Figs. [H and [8] are 
similar if the experimental limits and differences of cutoff parameters are neglected. In fact, 
we have checked that the patterns of the scatter plots in Figs. [H and [HI become the same if 
we draw all the sample points, including those excluded by the experimental limits as long 
as we take same cutoff parameters. 

Figures. [9] and [TOl show the asymmetries of the type II decay r+ — )■ /i^'^e+e^. As given in 
Sec. nil C| we define the seven asymmetries. We can see that the allowed ranges of the parity 
asymmetries Az, Ax, Azfb and Axfb are almost independent of the branching ratio. 
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FIG. 8: Angular asymmetries of r"*- — ?• fi^fj,~^fj,~ as functions of the branching ratio for the same 
parameter set as in Fig. [3j The correlation between Az and Ax is also shown in (c). Ay is 
magnified by 100 in (d). The vertical solid lines in (a), (b) and (d) are the experimental upper 
limit of Br(r^ — t- ^u^^^/i^). The color code is the same as in Fig. [3l 
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FIG. 9: Angular asymmetries of r"*" — )• fi~^e^e~ as functions of the branching ratio for the same 
parameter set as in Fig. [3j The correlation between Az and Ax is also shown in (c). Ay is 
magnified by 100 in (d). The vertical solid lines in (a), (b), and (d) is the experimental upper limit 
of Br(T+ — )■ fi'^e^e~). The horizontal dashed line in (a) is the value in A^!^'^ — t- limit. The color 
code is the same as in Fig. [3l 
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FIG. 10: Forward-backward asymmetry and forward-backward type angular asymmetries of — )■ 
/i"'"e^e~ as functions of the branching ratio for the same parameter set as in Fig. [3j Ayfb is 
magnified by 100 in (d). The vertical solid line in each plot is the experimental upper limit of 
Br(T+ —7- /i+e^e^). The horizontal dashed line in (c) is the value in A^"^ — t- limit. The color 
code is the same as in Fig. [3l 
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Taking the A^"^ — )• limit is informative to understand the behavior of the asymmetries, 
because the contribution of the dipole term in the LHT is relatively small. In this limit, 
values of Az and Axfb are constants: 

^ ' ALHT^O 105 ^ ' 

Deviations from these values seen in Figs. [9]^a) and fTOl c) are identified as effects of the dipole 
term. This effect is larger in Az because the function in Eq. fl62ap becomes logarithmically 
large for 5 — ?■ [see Eq. flA30eP ]. while the function in Eq. flA29cP does not have such 
an enhancement. On the other hand, A^^'^^"'" and A^^]^"'" depend on the relative magnitude 
of 5'^'^;'^^"'" and 5'^^^'"^"'", so that these asymmetries depend on input parameters even if the 
dipole term is neglected. That is why A-^ distributes within a wider range compared 
with A^''"^"'". 

Since the r+ — )■ yU+e^e^ decay in the LHT is described by only three Wilson coefficients, 
we can derive the following proportionality relation for 5 — )■ 0: 

Atx 11 

Ax : Afb : Azfb = • "4 • ^^^^ 

This relation holds for any values of 5'^!/'"^"'"; d^Lr'^^ and A^'^. Although ApB is not a parity 
asymmetry, it is related to the parity asymmetries Ax and Azfb because of the restricted 
chirality structure of the LHT. We can see that Eq. (!89|) holds in a good approximation in 
Figs. Mh),M,a) andHni^b). 

More relations characteristic to the LHT are discussed in Appendix O When the LEV 
processes are precisely measured, these relations might be useful to determine whether or 
not the new physics is the LHT. 

The time-reversal asymmetries Ay and Ayfb are very small for Br(r+ — )■ /z+e+e^) > 

10-13. 



2. Semileptonic decay modes 

The correlations between branching ratios of r"*" — t- /i+P (P = 7r°, rj, rj') and r+ — >■ 
are shown in Fig. [11] These branching ratios are roughly of the same order of 

I LHT 

magnitude for each set of input parameters, because relevant Wilson coefficients g£i ^ and 
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FIG. 11: Correlations between Br(r"'" — t- fJ-'^P) and Br(T"'- — > fi'^fi^fi~) for P = vr", r/ and r/'. 
Input parameters are the same as those in Fig. [3l The vertical and the horizontal solid lines are 
the experimental upper limits. The color code is the same as in Fig. [31 




9Lir{q) (HS]), dM]), (175]) and ( 171]) are similar in magnitude and the effect of the dipole 



term (in r+ — )■ on the branching ratio is small. 

In Fig. [12] we show the branching ratios of r+ — fi^V {V = p°, 0) decays. The 
qualitative behavior of the correlations is the same as in the r'^ — t- fi~^P case. Polarization 
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FIG. 12: Correlations between 'Bi{t~^ — t- /i"*"!^) and Br(r''~ — t- fi~^fi~^fi~) for 1/ = p'', cj and 
Input parameters are the same as those in Fig. [3l The vertical and the horizontal solid lines are 
the experimental upper limits. The color code is the same as in Fig. [31 

asymmetries are shown in Fig. [T31 The contribution of the dipole term in r+ — )■ fi~^V 
affects the polarization asymmetry though the effect on the branching ratio is smalL In the 
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FIG. 13: Correlations between the polarization asymmetries and the branching ratios for r+ — t- 
fi^V {V = p^, CO and 0) decay modes. Input parameters and the color code are the same as those 
in Fig. [3l The vertical solid lines are the experimental upper limits of the branching ratios. The 
dashed lines show asymmetries in A^"^ — )• limit. 



A^"^ — )■ limit, the asymmetry is determined as 



1ml — 2"^y 



4 



0.22 for V = p°, u, 
0.10 for V = 6. 



(90) 



2 m2 + 2ml 

The deviations from these values seen in Fig. [13] can be understood as effects of interferences 
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between the dipole and the four-Fermi terms. In particular, in the parameter region where 
the branching ratios are larger than 10~^°, contributions of the four-Fermi terms dominate 
the decay amplitude, so that values of Ay are not much different from those given in fl90p . 

As in the case of trilepton modes, we can derive several relations among various branching 
ratios asymmetries of semileptonic modes, which is discussed in Appendix O 

C. Correlations among — > e, r — t- ^, and r — > e transitions 

There are three classes of processes which change the lepton flavor by one: yU — )■ e, 
r — 7- /i, and r — )■ e transitions. Here, we discuss correlations among different lepton flavor 
transitions. 

We show correlations among branching ratios of /i+ — e^7, r+ — )■ /i^7, and r+ — )■ e"'"7 
in Fig. O As can be seen in Fig. [111(a) and (b), there are no direct correlations between 
the branching ratios of /i"*" — )■ e"'"7 and r LFV decays, because we take the T-odd lepton 
masses and parameters in the mixing matrix Vne as free parameters and vary them 
independently. On the other hand, in Fig. [TWc). we notice that a parameter region where 
both Br(r+ — )■ fi'^'j) and Br(r+ — t- e+7) are larger than 10~^ is not allowed. When both 
T ^ fi and r — )■ e transition amplitudes are large, the corresponding /i — )■ e amplitudes 
also become large so that the branching ratios of the /i — )■ e transition processes exceed 
the experimental bounds. In fact, we have checked that the branching ratios Br(^i — >■ 
£27)/Br(£i — )■ £21^1^) [{ii, £2) = ifJ', e), (r, /i) and (r, e)] distribute within the same range 
(< 10^^) when we ignore experimental constraints on yU — )■ e processes. 

The correlations among the type I leptonic three-body decay branching ratios are shown 
in Fig. [151 The behavior of the correlations is the same as the case of Fig. [TH 

V. CONCLUSION 

We have calculated branching ratios and angular and forward-backward asymmetries of 
fi and r-LFV processes in the LHT. We have obtained the following results: 

• The branching ratios of three /i-LFV processes, — e+7, fi^ —?■ e+e+e^ and the /i — e 
conversion rates, can be close to the present experimental bounds. There is a rather 
strong correlation between the branching ratios of /i"*" — )■ e"'"7 and fi~^ — )■ e~^e~^e~ . This 
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FIG. 14: Correlations among branching ratios of /x"*" — )• e"'"7, — )• /i^7 and r"*" — )• e"'"7. Input 

parameters and the color code are the same as those in Fig. [3l The horizontal and the vertical 
lines are experimental upper bounds. 

is in contrast to the case of fi^ — t- e+7 and theu — e conversion where no correlation 



is observed. These features are noted in Ref. 



12| 



The parity asymmetry of e"'"7 is -1/2. For /i+ — )■ e~^e~^e , the parity asymmetries 

can be large, whereas the time-reversal asymmetry is within 10%. 
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FIG. 15: Correlations among branching ratios of fi~^ — )• e~^e~^e~ , — )• ^'^ ji'^ ji~ and r+ — )• e^e'''e~. 
Input parameters and the color code are the same as those in Fig. [3l The horizontal and the vertical 
lines are experimental upper bounds. 

• We have calculated the /i — e conversion rates for various muonic atoms: Al, Ti, Au 
and Pb. In most of the parameter space, ratios of the conversion rates are found within 
1 order of magnitude. In some cases, however, the conversion rates for Ti, Au and Pb 
can be close to the experimental bounds even if the rate for Al is suppressed. 
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• The maximal values of the branching ratios for various r-LFV processes are 10~^-10~^ 
except for 0(10^"'^°) for r+ — t- yU+w and r+ — )► e^cu, and 0(10^^'^) for r+ — t- fi^fi^e~ 
and r+ — )■ e+e+/x^. 

• The parity asymmetries of r"*" — t- /i+7 and r"^ — > e+7 are -1/2 and r"*" — )■ yU+P and 
r+ — !■ e+P are 1/2 reflecting the chirality structure. For r+ — )■ /i+y and r+ — > e~^V, 
if branching ratios are larger than 10"^*', the asymmetries are about 0.3, 0.15, and 0.1 
for u, and 0, respectively. 

• There are sizable parity asymmetries in r+ — )■ fi^fi^fi^ , r+ — t- e+e+e~, r+ — )► yU+e^e~, 
and r+ — )■ e'^fi'^fi". For r+ /i+e+e^ and r+ — )■ e^^^fi^, forward-backward asym- 
metries can be defined, and there is a relation among asymmetries as Eq. flS^ . The 
time-reversal asymmetries are found to be very suppressed. 

The search for LFV in fi and r decays has a complementary role to the new particle search 
at the LHC experiment to explore the TeV-scale physics. We have seen that branching ratios 
of these processes in the LHT can be within the reach of ongoing and planned experiments 
such as MEG, COMET and Mu2e for /i-LFV processes and r rare decay searches at future 
B factories and the LHC. We have found that various asymmetries defined with the help of 
/i and T polarizations reflect the characteristic chirality structure of the LHT. Experimental 
searches for LFV processes using polarization of initial leptons are important to identify this 
model among various candidates of TeV-scale physics models. 
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Appendix A: General formulae of branching ratios and asymmetries 



We present general formulae of the branching ratios and asymmetries in the LFV processes 
based on the general low energy Lagrangians. Some basic formulae can be found in Refs. [l^ 



ml4i|. 



1. —7- T+ —7- and r+ — t- e^7 

The Lagrangian of the radiative two-body decay is 
AGf 

C'y = -^[mrARfRa'^yLF^,, + mrALfLa''yRF^, + H.c.]. (Al) 
The differential decay width of — t- is 

dBr(rJ ^1.--,) ^ r^9>^^^A,\^ + |.4«p + (|.4,P - |^hP)co=9). (A2) 

a cos U TT 

The branching ratio and the asymmetry defined in Eq. fl4Up are 

Br(r+ ^ /.+7) =r.^^(|A^P + I^^H, (A3) 

TT 

A(^"->/^"7) 4 I1"[:!^[ - (A4) 



2. /i"*" —7- e"'"e"'"e , r"*" — t- fi~^fi~^fi and r+ — t- e"'"e^e 

The Lagrangian of the type I leptonic three-body decay is 
Ci = — [mrARfRa^" fiLFf,^ + mrALfL(y^''^RF^y 

+ 9Rri'^Rl^f^R)if^Rlt^f^R) + 9Lli^Ll^f^L){fiLl^f^L) 

+ 9Ri{^Rl^m){fiLl,,f^L) + ^L(TL7>L)(/iiJ7M/^i?) + H.c] . (A5) 
The differential partial decay width of r"*" — fi~^^'^fi~ can be written as 



TV 



{MQ{xb,Xc) + Ml{xb,Xc) COS 6 



dxbdxcdcj) dcosO 1287r*^ 

-|- Mx{xb, Xc) sin 6' cos -|- My(a;b, Xc) sin 6* sin 0), (A6) 
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where 



Mo{xb,Xc) =(C]ji + Cl^)ai{xb,Xc) + (C]^2 + Cl2)a2{xb, Xc) 
+ iC]i^ + Cl3)a3{xb,Xc) 

+ (Cji + Cj2)a4(xb, Xc) + (Cj3 + Cj^)a5(xb, Xc), 

Mz{Xb,Xc) =iC\i - Cl^)bi{xb, Xc) + (C]j2 - Cl2)b2{Xb,Xc) 

+ {CR3-Cls)a3{xb,Xc) 

- (Cji - C]2)a4(Xb, Xc) + (Cj3 - Cj4)a5(x6, Xc), 

Mj^{Xb,Xc) ={Cri - Cl,)ci{Xb,Xc) + (C^2 - Cl2)c2{Xb,Xc) 

+ - Cj2)c3{Xb, Xc) + {Cjs - Cj^)c4{Xb, Xc), 

MyiXb, Xc) ^Cj^Cs{Xb, Xc) + C]QC4{Xb, Xc), 



The functions are defined as follows: 

fl - Xb){2xl - 2xb + !) + (!- Xc){2xl - 2xc + 1) 



ai{xb,x, 

a2(xb,x, 
a3{xb,x, 
a4{xb,Xc 
a5(xb,x^ 



(1 - Xb){l - Xc) 
= 2{xb{l - Xb) + Xc{l - Xc)), 
= 8(2 -Xa- Xb){Xb + Xc-1), 

= 32{xb + Xc- 1), 
= 8(2 -Xb- Xc), 



(A7a) 



(A7b) 

(A7c) 
(A7d) 



^Ri, -- defined as 








C],,^\eAnf, 




^\eAL\\ 


(A8a) 


Cr2 = Ifl'ral ) 




— \9Lr\ 1 


(A8b) 


_ \9rs\ i t |2 
^R3 ~ IQ \9Rr\ 1 




1 I |2 

|5'L6-| , I I |2 

- -J^g + Ifl'wl > 


(A8c) 


C]^ = MeARg'li], 




= MeALo'^X 


(A8d) 


C}3 = ReM«^i;], 




^Re[eALg%], 


(A8e) 


C)^ = lm[eARg^li + eALgtrl 




= lm[eARg^l^ + eALg%]. 


(A8f) 



(A9a) 

(A9b) 
(A9c) 
(A9d) 
(A9e) 



{xb,Xc) = 8(^ 



l-2xb{l-Xb) ^ 1 



2Xc{l-Xc) ^{Xb + Xc-l) 



Xr 



I- Xb 



Xj) Xq 



). 



b2{Xb,Xc 



{xb + Xc){xl + xl- 3{xb + Xc) + 6) - 4 



Xi) Xc 



(AlOa) 
(AlOb) 
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ci{xb,Xc) = -32 J- -, Alia 

2-Xb-Xc y {1 - Xb){l - Xc) 



C2{xb,Xc) = -4:-^ — ^^a/(1 - Xb){l - Xc){xb + Xc- 1), (Allb) 
C3{xb,Xc) = -16{xb - Xc){xb + Xc - ^)\ TT^^-YT, — (Allc) 

y (1 - Xb){i - Xc) 

C4{xb,Xc) = -S{xb - Xc)(2 - Xfc - Xc)\ —^^-^^^ — (Alld) 

y (1 - Xb){i - Xc) 

where Xb and Xc are defined in Eq. (H7|) . The branching ratio and angular asymmetries 
defined as Eqs. f l51ap - (l51cp are 

Br(r+ ^ =Br(r+ ^ z/,e+z/e)S^(5), (A12) 



+ (C'n + C'jM^iS) + (C^s + C'jM^iS), (A13) 



+ iC],, - CiMsiS) - [C], - C],)A,{5) 

+ {C],-C],)A,{5)), (Al4a) 

= ^^^{-(Cri - Cli)Ci(^) - (Cm - Cl,)C2{5) 

- (Cji - C'j,)C,{6) + {C'j, - C],)C,{6)), (A14b) 
A'viS) = ^^^i-Cj^CsiS) + C'j,C,{6)). (Al4c) 

Br, A\, A\ extract the components Mq, M\ and My in Eq. flA6p . respectively. 

Since the signs of cos^, sin 6^ cos and sin 6* sin are equal to s ■ pa, s ■ {{pa x pb) x pa) and 
s ■ (pa X p;,), respectively (see Fig. [2]), and A]^ are parity odd asymmetries and Ay is 
a time-reversal asymmetry. The following functions are introduced in the above formulae. 
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The cutoff parameter, S, is defined in Sec JIIII for each of the processes. 

Ai{6) = -16(1 -6){2-6 + 26^) In (^j^) ^'^^^^^ " ^'^ + 

^2(5) = {l + 26 - 26'){l-26f, 
A,i5) = 2{l + 25){l-25f, 
A^{6) = 16(1 -25)^ 
^5(5) = 8(1 + (5)(1- 25)2, 



Bi{5) = - 16(2 + 215 + 35^ - 25^) \n{25) + 16(1 ~5){2-5 + 2S^) ln(2(l 

o 

- -(l-2(5)(49 + 685 + 452), 
S2((5) =^(1 - 2(5)(1 + 86- 386^ - 126^) - 166^ \n{2S). 



Ci(5) = — (4 + 96 + 6^)VT^^ - 48v^(3 + 66- 6^) arccos ( - 
5 V 1 



36 -r 



3846 arccos 



1-6 



4 



C2(6) = (8 + 86 - 936^ - 2256^)Vl-26 

105 

3 /' 36 1 
— 26^(1 — 65 — 3(5^) arccos I — ) — 165"^ arccos 



1-6 \l-6 



Cs(6) = — 71^25(48 - 575 - 685^ + 855=^) - 12^5(1 - 6)'^ arccos ( ^^-^ ) 
35 \ 1 — 6 J 

Ci{6) =—71^25(64 - 415 + 265^ - 855^) 
35 

- 6^5(1 -6-6^ + 6^) arccos (yrj) • 
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3. r"^ — >■ /x+e'^e and r"*" ijl 

The Lagrangian of type II leptonic decay is 

+ 9Rsi^Rf^L){eReL) + gLs{fLf^R){eLeR) 
+ 9Rt{^ReL){eRHL) + 9Lti^LeR){eLfiR) 

+ 9RriT'Rl^l^R)i^Rln(^R) + gLli^Ll^l^L)ieLlt,eL) 

+ 9Ri{^Rl^f^R){eL7t^eL) + 9Lr{rLl''f^L){eRlt,eii) 

+ 9iL{rR-f''eR){eL^^fXL) + 9ii{rLl''eL){eR^^lJiR) + H.c] . (A18) 
The differential partial decay width of r"*" ii'^e^e~ can be given as 



d^Bi{T+ /i+e+e") Glnv> 



[Ml{Xb, X,) + (Xfe, Xe) COS^ 



dxbdxc d(f) dcos 9 12 Stt*^ 

+ Mx{xb, Xc) sin ^ cos + Myix^, Xc) sin 9 sin 0), (A19) 



where 



M'^{Xb, Xc) ={C'^i + C'ii)di{Xb, Xc) + + CL2)d2{Xb, Xc) 

+ + C]^s)ds(x,, Xc) + (C^, + C'j^dd4(x,, Xc) 

+ (Cji + C'jMx,, X,) + {C'J, + C'J,)d,{x,, Xc) 

+ {C'I, + C%)dr{x,,xc), (A20a) 

M'iix,, Xc) HC'^, - Cl\)ei(x6, Xc) + - C'^2)d2iXb, Xc) 

+ (^«3 - C'^3)e2{xt, Xc) + {C^, - C'^desixt, Xc) 

+ {Cf, - Cf,)e,{x,, Xc) + {Cf, - C'J,)e,{x,, Xc) 

+ (C'J,-C'jMx,,Xc), (A20b) 

M'^{x,,xc) ={C'i, - C]^,)Mx,,Xc) + {C^, - C]^,)f2{x,,Xc) 

+ {C'^4-C]^dMxb,Xc) 

+ {C'k - C'J^)h{xb, Xc) + (C^a - C'j,)h{x,, Xc) 
+ iCl-C%)Ux,,Xc), (A20c) 
M^{xb, Xc) =Ckgi{xb, Xc) + Cfs92{xb, Xc) + Cfgfeixb, Xc). (A20d) 
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C]ji ... are defined as 



Cl^^\eAR\\ Cl^^\eAj\\ (A21a 

Cl,= ^^^\9n.\\ Cl,J^^\g,.\\ (A21b; 

Cl,-^^^\9m\\ Cl,J^^\g,r\\ (A21c 

Cl^^\9Rr\\ Cl^^\gu\\ (A21d: 

C\ = ReM«^2,], = ReMi^^,], (A21e 

= ^^{eAnglA. C'j, = MeAL9R,\, (A21f; 

= Reb^,^;^,], = Re^^.^^t], (A21g: 

= Im[eA^y2/ + (^^LgRr] , C-^g = MeAngl^ + e^^y^J , (A21h; 

The functions, di{xi),Xc), ei{xb,Xc), fi{xb,Xc) and gi{xi„Xc), are defined as 

Xfe + OJc - 1 

«i2(a;b, Xc) =2(2 -Xb- Xc){xb + Xc- 1), (A22b^ 

dsixb, Xc) =2x6(1 - Xb), (A22c' 

d4^{xb,Xc) —2xc{l — Xc), (A22d 

d5{xb, Xc) =8(1 - Xc), (A22e' 

d6{xb,Xc)^8{l-Xb), (A22f 

d7(x6, Xc) = - (1 - Xb){Xb + Xc-l), (A22g^ 



^l('^6) Xc 

e2{xb,Xc 
e3{xb,Xc 

^4{,Xb, Xc 
^bi^Xb, Xc 



Xb + Xc-1 



Xb{l - Xc) + Xc{l - Xb) -2 



(1 - Xb)' + (1 - Xc)' 
2 Xb Xc 



2 Xc 

V 2 — Xh — Xn/ 

= -8i 



2 

,(1 - a;c)(xfe - Xc) 



=8 



2 Xb Xc 

(1 - X{,)(Xfe - Xc) 

2 Xj, Xc 



(A23a) 
(A23b) 
(A23c) 
(A23d) 
(A23e) 



48 



fi{x,,x,) =16 /^ "2^^^ (A24a) 

2 - Xb - Xc\j Xb + Xc- I 

f2ixb,Xc) =4:-^ — V(l - Xb){l - Xc)ixb + Xc- 1), (A24b) 
h {xb, Xc) = - 4 ^""^' y/{l-Xb){l-Xc){xb + Xc-l), ( A24c) 

Z Xjj Xc 

+ (A24d) 

I — Xb — Xc y Xfe + Xc — 1 

J = - 3 (1 - + X,) /a - 

/ — Xfe — Xc y Xb + Xc — 1 

/6(Xb, Xc) = - V(l-X;,)(l-Xc)(Xb + Xc-l), (A24f) 



^i(x,,xc) = - 24(1 - xc)\ r. """^^^ (A25a) 

y (xb + Xc - 1) 

<72(x,,Xc) =24(1 - Xb)^!^^^^^. (A25b) 
The branching ratio and seven asymmetries defined in Eqs. f l51al) - fl51cl) and flMl) - fl65cp are 

Br(r+ ^ /x+e+e-)(5) =Br(r+ ^ z/,e+z/e)5"(5), (A26) 



+ (^^"3 + C^ia + Cl, + ClijD3(5) + {C'X + + C^J3 + C'k)m). (A27) 
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+ {CI,-C'l,-\{C'J,-C'JS)D,{5) 

- (Cr3 - C'is + C^m - CLi)D&{5) 

- \{C'k - C'J, + C'k - C'h)D,{5)^ , (A28a) 

+ (C'A - C'J^ - C'J, + C'J,)E,iS)^ , (A28b) 

A'^iS) = (-(Cy, - C'Js)Es{5) - lc'l,E^{6)^ , (A28c) 

+ \iC\ + C\ - Cl - Cf,)D,{6)^ , (A29a) 
A'iMS) = (-^ (c'i, - C'i, - C'^, + C'^, - ^-{C'J, - C'J,))d,{S) 

- liC'A - C'J, - C'h + C'J,)D,{5)^ , (A29b) 

+ liC'J, - C'J, + - C'h)E,{5)^ , (A29c) 

AvFBiS) = (-^(Cy, + C'J,)E,(5)^ , (A29d) 
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with 



16 

Di{5) = —(-(1 - 5){S -S-5^)- 6\nS), (A30a) 

02(6) = (1-6)^1 + 35), (A30b) 

D3{6) = {l-6)%1 + S), (A30c) 

D^{5) = 8(1 - 5)^ (ASOd) 
16 

^5(5) = —(-(1 - 5)(10 - 55 + 5^) - 61n5), (A30e) 
3 

L'6(5) = ^(1-5^(1-35), (A30f) 



EUd) = — (1 - 5^/2)3/g ^ + + 45^3/2 ^ ^^^2^ .^3^ X 

35 

£;2(5) = — (1 - 5^2)3(64 + 875^2 + 695 + 455=^/" + 155"), (A31b) 

Es{S) = 7^(1 - 5^/2)^(16 + 295^/2 + 205 + 56^^^), (A31c) 
35 

32 

£;4(5) = y-(l - 5^/')'(8 + 95^/2 + 35). (A31d) 



4. r+ ^ ju+^i+e and r+ ^ e+e+/i 

The effective Lagrangian of type III leptonic decay is 

An = -—^[gnJ {fRiiL){eR^L) + gLs\fLiiR){eLiiR) 

+ 9RrifRj^HR){eRj^HR) + giYifLl^HL)ieLl„I^L) 

+ QmirR^I^R) (eL^t^l^L) + gir {tl^'' I^l) {eR^i^iiR) + H.c] , (A32) 

The differential partial decay width of r"^ — )> /x^/x+e" is written as 

dxbdxcd(pdcose ^ 12871^^ o \ b, z \ b, o 

+ {xb, Xc) sin 9 cos sin ^ sin 0), (A33) 
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where 

M'^\x,, ={C'^l + C™)«2(X6, + {C'^l + C'^l)as{x,, xj, (A34a) 
M™(x,, Xe) =(C™ - Cli2')&2(x6, X,) + (C'^l - C]^l)as{x,, xj, (A34b) 

M™(x6,Xe) =(C]jI^ - Cll2')c2(Xb,Xe), (A34c) 

M^"(xfe,Xe) =0. (A34d) 
Integrating them in x^-Xc plane, branching ratio and parity asymmetries are given as 

Br(r ^ =Br(r+ ^ P.e+z.e)((C]ji^ + Clis')^^!^) + (C],^^ + C^l)As{6)), (A35) 

^ni/.N (Cfl" - Ci2)C2(^) . . oyN 

^^^'^ 2(C7],UC^i^2^M2(5) + ra + C^i^i)^(5)' ^ ^ 

=0, (A38) 

where coefficients Cj^^s are given by replacement of I with III in coefficients of the type I, 
Eqs. flA8ap - flA8f|) . In this type, time- reversal asymmetry. Ay, does not appear. 



5. Semileptonic decays of r 

The effective Lagrangian of semileptonic r decays is 

aGf r 



q=u,d,s 

+gRt(q)ijRlJ,L){qLqR) + 9Lt(q)ifLlJ,R)iqRqL) 

+9Rr(q){fR-i'' ^R){qR-i^qR) + gLi{q){fLl^HL){qLl,,qL) 

+9Rl{q){rRl^f^R){qLlt,qL) + 9Lr{q){fLYf^L){qRlt,qR) 

+9RT{q){fR(y'"' fiL)iqRcrfiuqL) + 9LTiq)ifLa'"' fXR){qL(y^,uqR) + H.c.) . (A39) 
The branching ratios and asymmetries for r+ — )■ fi'^P decays are 

Br(r+ ^ = f 1 - 4V d^^^l' + l^^^l')' (^40) 

Att \ mi / 



AAr- = i^^^;^!^. (A41) 
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Here, Grp and Glp for P = it and rj are given by 

^^{gRt{u) — gRt{d) — gRs{u) + gRs(d))j , 



2m. 



(A42a) 



Gltt ~ir7K {gRr{u) ~ gRr{d) " gRl{u) + gRl(a 



2^2 



+ \gLt{u) — gLt{d) — gLs{u) + gLs{d)) ) : 



(A42b) 



1 / /9 

Gr^ =- \^--^{gLr{u) + gLr(d) " " + f^igLr{s) " 



2^2; 



{gRt{u) - gRt{d) - gRs{u) + gRs{d)) + {gRt{s) - gRs(s)) ) , (A42c) 



rnqTiij- 



2m„mr 



.1(11 



Glt) —- y-^i.gRriu) + gRr{d) — gRl{u) — gRl{d)) + f^igRr{s) — gRl{s)) 

\gLt{u) — gLt{d) — gLs{u) + gisid)) + 7^ [gLt{s) — gLs{s), 



2mgm^ 



(A42d) 



2\/2m 

The expressions for rj' are obtained by the replacement — )■ r/' in the last two equations. To 
describe the decay constants for rj and r^', we used the formalism in Ref. 42[. We assume 
isospin symmetry so that m„ = = iriq in the following. The decay constants are 



-iUp^ = (0 

-^UP' = (0 

-^fsP^ = (0 
(0 

(0 



2mq 

ihs_ 
2m, 



V2_ 



V2 



VqiP)), 



V2 



\Vgip)), 



sj^s\risip)), 



(A43a) 

(A43b) 
(A43c) 
(A43d) 

(A43e) 



For T] and r]', the decay constants, /^o'^''' are h^o'^'^ are expressed by fq, fs and a mixing 
angle, (pri, as follows. 



fq fs 

J r) J rj 

fQ fs 

Jri' Jri'. 



COS (Pri — Sm (Pri 



Sm (pr, COS (pr, 



fq 
/. 



(A44) 



COS </)^ — sm 



sm {/)^ cos (/)^ 



hq 

/i. 



(A45) 
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where 



hq — fq{n^l cos^ (j)r, + 171^, sin^ (prj) — V2fs{m'^, — m^) cos 0^ sin 0^ (A46a) 
hs = fs{ml, cos^ 0^ + ml sin^ 0^) - -^fs{ml, - mj) cos 0^ sin 0^. 



(A46b) 



Similarly, for vector mesons, branching ratios and asymmetries are 
Br(T+ ^ /i+y) 



^2 _2 I o^,^ 



X ( Cav+ + <-y+ — 6Re[Liv+\ ] , (A47j 



, . + . ,^+.,^ ^1 C'Ay-(2m^ - ml)/ml + Cy_(m^ - 2m^)/(4m^) + Re[Cjy-] 

^ ^ 2 CAy+(2m2 + ml)/ml + Cy+(m2 + 2ml) /{Ami) - 3Re[C,y+] ' ^ ' 



with 



(A49a) 
(A49b) 

Civ+ — GravG*j^y + GlavG*jiy, Gjv- — GRAvG*i^y — GlavG*jiy. (A49c) 



Cav+ — 


\GLAvf + 




Cv+ — 


\Glv\^ + 


\Grv\^ 


Cav- = 






Cv- = 




\Grv\^ 
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The effective couplings Gray, Grv, Glav and G^y for y = a; and are 

^ _ Jp 9rt{u) - 9RT(d) , f rrir Qu-Qd ^^ ( 
GuAp - - Y ^ + U—^^eAn, (A50a) 

^ fj 9LT(u) - 9LT{d) ^ fTTlrQu-Qd . , . x 

G,Ap = + U--^-^^^ (A50b) 

f 

^RP = ITIk (9Rr(u) - 9Rr(d) + 9RI(u) " 9Rl(d)): (A50c) 

2v2mr 

^Lp = —7= {9LI{u) - 9Ll{d) + 9Lr{u) " 9Lr{d)), (A50d) 

2\/ ZlTLr 

GuA. = Jl 3mu)+_9RTid, ^ ^v^Qu±Qd^^^^ (^50^) 

G.^. = + iJ^'^^eA,^ (A50f) 

2 V2 V2 

^ n 'at + + + fi-fflCd)), (A50g) 

2v2mT- 

^Luj = -r^7^{9Ll{u) + 9Ll{d) + 5'Lr(«) + 9Lr{d)), (A50h) 

2v2mr 

Gra,/, = — ^9rt{s) + f4>—QseAR, (A50i) 
Gla,^ = — ^9lt(s) + /0— Q^e^L, (A50j) 

z 777^ 



Gr4> = -^^idRris) + gRl{s)): (A50k) 

GL^^^{9Llis)+9Lris)) (A501) 

The decay constants fy are defined as 

Tn./pe'^ = (o r ^ ^ (A51a) 

m./.e'^ = (O r^'^J^^^'^ |a;(p)), (A51b) 

7770/06'^ = (O|s7'^s|0(p)), (A51c) 

where e'^ is the polarization vector of the vector mesons, fy are also defined with the same 
flavor combinations, 

m<^'''<l\y{p)) = -^/y (P'^e^ - P""^")- (A52) 
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Pseudoscalar 


Vector 


u 


130 MeV 




221 MeV 


fq/ fir 


1.07 


fuj 


196 MeV 


fs/ fiv 


1.34 


f(f> 


228 MeV 


4>V 


39.3° 





TABLE III: Decay constants of pseudoscalar and vector mesons. 

fv can be extracted from partial decay widths for e~^e~ decay modes, 

J p 
f 

J Ul 

f 2 _ 

The values of the decay constants are listed in Table IIIII 



p2 3T7lprp_^e+e- 

27?7T.(jr(^_^e+e- 

27ra2 



(A53a) 
(A53b) 
(A53c) 



6. fj, — e conversion 



The effective Lagrangian for /i-e conversion is obtained by the replacement of r — )■ /i and 
— )■ e in Eq. ( ]A39I) . According to Ref . j22| the rates of coherent /i — e conversions are 

2 



R{H~A -> e~A) 



2Gl 



l^capt 

+ {L^R) 



-A,D + 2 y: + -afM""' + CAT + ~9iV)vn 

N=p,n 



(A54) 



where Wcapt is the muon capture rate of each atom and the p and n mean proton and neutron, 



respectively. 



in the above equation are written in terms of the Wilson coefficients 



in the effective Lagrangian: 

9lI + 9Lt = Y '^S'^^ {dLsiq) + 9Lt{q)), 
Q 

9ls + ~9U = Yl ^5'"^ {9Ls{q) + 9Lt(q)): 



~ (p) ~ (p) 

9 LI + 9lt = '29U(u) + '2gLr{u) + gLl{d) + fl'Lr(d), 
9^1 + 9^} = 9LI{u) + 9Lr{u) + 2to(d) + '29Lr{d), 



(A55a) 

(A55b) 

(A55c) 
(A55d) 
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where G^'^-* are given in Ref. 22|: 





_ Q(d,n) 


= 5.1, 


(A56a) 




- '^S 


= 4.3, 


(A56b) 




= Gg'"^^ 


= 2.5. 


(A56c) 



D, S^P\ V^P^ and K^"), in Eq. (IA54D are overlap integrals defined in Ref. j22j. They 
depend on nuclides. In this paper, we calculate the conversion rates for Al, Ti, Au and Pb 
using the numerical values listed in Table llVi 



Nucleus D F^^) Wcapt 

llA\ 0.0357 0.0153 0.0159 0.0163 0.0169 0.7054 
||Ti 0.0870 0.0371 0.0399 0.0462 0.0495 2.59 
ifjAu 0.167 0.0523 0.0859 0.0610 0.108 13.07 
2^fPb 0.163 0.0493 0.0845 0.0686 0.120 13.45 

TABLE IV: Values of overlap integrals and capture rates. They are taken from Table II for Al 
and Au, IV for Ti, VI for Pb and VIII in Ref. [23] • The integrals are in units of m^^'^ and rates are 
m units of 10 V sec. 



Appendix B: Functions used in the Wilson coefficients 

The following functions are used to describe the Wilson coefficients in Sec lIIII 

^cm{x) =— ^ -1 + 5a; + 2x^ + Inx , (Bla) 

4(1 — xy \ 1 — X / 

Nnm{x) .3 (2 + ^x-x^ + In , (Bib) 

4(1 — xy \ \ — X J 

P^{x) = - 2 (^Ncc{x) + ]^Nmc{x) + ^Nmc{x)^ , (Blc) 

2x \ 

12 + a;-7x2 + (12-10a; + x2) Ina; , (Bid) 

1 — X J 

2x \ 2 

18 - llx-x^ + (15 - 16x + 4x^) Inx - -Inx, (Ble) 

X X y 3 

6-x+(2 + 3x)-i^ ) , (Blf) 



Ncc 


[x) 




X 




12 


[1- 


xY 


Nnc 


[x) 




X 




12 


[1- 


xf 


Pz 


[x) 




X 


-{ 




1 - 


X \ 



1 — X 
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B{e){x,y) =- {Bc{e)ix,y) + BN{x,y) + BNx{x,y)) , (B2a) 
B(^u){x,y) =i {Bc{u){x,y) + BN{x,y) - BNx{x,y)) , (B2b) 
B(d}{x,y) =^ {Bc(d){x,y) + BN{x,y) + B,yx{x,y)) , (B2c) 



Bc{u){x,y) ={16 + xy)Bxix,y,l) - 8xyBxc{x,y), (B3a) 
Bc{d,e)ix,y) = - (4: + xy)Bxix,y,l) + 8xyBxcix,y), (B3b) 
Bxix, y) =3 {Bx{x, y, 1) + ^Bxiix', y', 1) j , (B3c) 

BNx{x,y) =—Bx{x,y,r]) (B3d) 
5 

rj'^lnr] x'^lnx 

Bx{x,y,V) = 



{1 - ri){r] - x){ri - y) {1 - x){'r] - x){x - y) 



{1 - y){T] - y){y - x)' 

1 /xlnx ylny 

(l-x)(l -y) ' x-y v(i-a;)2 (i-y)^ 



Bxcix, y) =7^-37^ + ::— ( tt-^ - > (B3f ) 



where, 



(B3e) 



x = ^, x' = ^ = ^, ^ = ^^. (B4) 

Appendix C: Consistency test of the LHT using LFV branching ratios and asym- 
metries 

We introduce useful formulae to parametrize the various branching ratios and asymme- 
tries in the LHT. Since the observables are expressed by a small number of parameters, we 
can obtain implicit relations among these quantities, which can be used for a consistency 
check of the model. 

The nonzero elements of the Wilson coefficients are only of three types: 9^'^ 
and gYi^ ■ Here we consider the r — )■ yU transition. Corresponding formulae for the r — )■ e 
transition are also derived. According to Eqs. (145|1 . ( 146|) . fl58|) . ( I73|) and flTijl . we can derive 
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the following parametrization: 



^I,LHT g^'^^^ 

g^LHT = ~ ^' g^LHT = - ^ + ^ , (Cla) 



II.LHT II.LHT 

^Lr __ ^ id 

e^LHT ^' g^LHT 

„LHT «LHT 



^-z + A", (Clb) 



z is a complex-valued function associated with photon penguin diagrams and a part of Z 
penguin diagrams. 

In leptonic decays, the following formulae are obtained using the above parameters: 

a Br(T+ ii+-f) \ 3 J ' ' 

+ 8Re[-3^ + 2A^], (C2a) 

-4(41n(5) + f)+l|.r-|-. + Af 

-4Re[-z + 2A^], (C2b) 

|(^-21 + lH^-Re[-z + 3A^]), (C2c) 

a Br(r+ — )■ /x+7) 35 

967rBr(T+ u+e+e-) / 4\ , ,2 , ,ni2 ^ r a m 

(C3a) 



a Br(T+ — >■ /i+7) 


967rBA^(r+ 


fi^fi^fi ) 


a Br(r^ 


- ^ /i+7) 


96n BAx(r+ 


— )■ fi^fi~^iJ,~) 


a Br(r^ 


- -> /i+7) 


967rBAy(r+ 





q; Br(T 


■+ ^ //+7) 


967r BAz(t+ 




a Br(T^ 


- ^ /x+7) 


967r BAx(t+ 





a Br(T+ //+7) 35 
967rBAy(T+ //+e+e-) _A8tt 
~a Br(T+ //+7) ~"35" 



(C3b) 

(^- 1^1' + \-z + A"f + 16Re[A"]^ , (C3c) 
Im[A"], (C3d) 
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— ^^#^^4^ =-(- \A' + I- + + 16Re[A"]), 



a Br(r+ — )■ /i+7) 
967rBA^i.B(r+ ^ ^u+e+e-) 

a Br(r+ — t- /U+7) 
967tBAxfb{t^ /x+e+e^) 
a Br(r+ — )■ 
967rBAyi.ij(r+ ^ /x+e+e") 
a 



(C4a) 
t 

- ^ (- l-^l' + 1-^ + A"| + 16Re[A"]) , (C4b) 

o 

(56+ + |-2 + A"|' + 16Re[-2z + A"]) , (C4c) 

(C4d) 
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fi4 

Im[-2;z + A" 
35 ^ 



Br(r+ yU+7) 

where BA^ (a = Z, X, Y, FB, ZFB, XFB and YFB) denotes asymmetry, Aa multiplied 
by the branching ratio. We take the limit 5 — )■ except for logarithmically divergent terms 
in Ai{6), Bi{6), Di{6) and -05(5). The proportional relation ( l89l) is found from Eqs. (EHcj), 
(1C4al) and (lC4bl) . Furthermore, ten independent observables can be expressed by three 
complex functions. 

Similarly, the following formulae are derived for semileptonic decays, 

16 Br(r+ ^ fi+n^) irf^ 



a Br(r+ 


^ /i+7) 


ml 


1 Br(r+ 


^ 




2a Br(r+ 


^/i+7) 


ml 


1 BA(r+ 






a Br(r+ 


^/i+7) 


ml 



\Au-Ad\\ 

(1 - 
Xp{l + 2xp) 

l + 2xp 



(C5a) 



+ 



32 

(1 - a; 



2z + /\., 

\2 



A, 



12 



1 + 2x. 



(C5b) 



+ 



1 - 2x 



32 



2z + A„ - Ad + 



2Xr 



(C5c) 
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Br(r+ ^ ^y)) _ vr/,^ 



a Br(r+ 




ml 


1 Br(r+ 


1 \ 


_ ^fl 


2a Br(r+ 






1 Br(r+ 




el 


2aBr(r+ 






1 BA(r+ 






a Br(r+ 






1 BA(r+ 






a Br(r+ 


^/i+7) 





1 



Jit 



+ 



l + 2x. 



32 



2 4 
-;z + A„ + Ad 

3 l + 2x.. 



'1-x ^2^2jl-x^ 
\9x^{l + 2x^ 



+ 



1 + 2x<^ 
16 



2 . 4 

3 ' l + 2x^ 



1 - 



1 



+ 



1 - 2x 



(l-x^)' 



32 
2fl 



9x^(1 — 2a;^ 



+ 



l-2a;„ 



16 



- A 4 

3^^ ' 3(1 -2x<^) 



(C6a) 



(C6b) 



(C6c) 



(C6d) 



(C6e) 



where Xa = m'^/m'^. From the above expressions, we can see that various consistency tests 
are possible with complex functions z, A„ — A^, A^ + A^ and A^. 
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